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Abstract 

Electronic and transport properties of Graphene, a one-atom thick crys- 
talhne material, are sensitive to the presence of atoms adsorbed on its surface. 
An ensemble of randomly positioned adatoms, each serving as a scattering cen- 
ter, leads to the Bolzmann-Drude diffusion of charge determining the resistivity 
of the material. An important question, however, is whether the distribution 
of adatoms is always genuinely random. In this Article we demonstrate that 
a dilute adatoms on graphene may have a tendency towards a spatially corre- 
lated state with a hidden Kekule mosaic order. This effect emerges from the 
interaction between the adatoms mediated by the Friedel oscillations of the 
electron density in graphene. The onset of the ordered state, as the system is 
cooled below the critical temperature, is accompanied by the opening of a gap 
in the electronic spectrum of the material, dramatically changing its transport 
properties. 
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When can an apparently random system be considered ordered? Or can an ap- 
parently random ensemble of impurities in a system be correlated enough to force 
the reconstruction of the electronic band structure in a material? In this Article 
we predict that a dilute ensemble of adatoms sprinkled randomly over a graphene 
monolayer [Tl, l2| can establish long-range correlations between their positions, despite 
the fact that they may be many graphene unit-cell lengths apart. This correlation is 
strong enough that at a transition temperature it will induce an energy gap in the 
electronic spectrum despite the fact that, in the "ordered" state, the distribution of 
adatoms does not show any crystalline structure. It rather resembles the ferromag- 
netically ordered state of spins of magnetic ions in dilute magnetic semiconductors [3j. 
The physical mechanism behind this phenomenon is the electron-mediated interaction 
between the adsorbents, which prompts their partial ordering into a configurations 
associated with a hexagonal superlattice with a unit cell three times bigger than 
that of graphene. Since the density, p of adsorbents is low, they occupy a small 
randomly chosen fraction of the equivalent positions on the superlattice. This order- 
ing folds the Brillouin zone and thus opens a spectral gap for low-energy electrons. 
This phenomenon suggests a novel route towards engineering the band structure and 
controlling transport in graphene-based devices. 

Graphene [H E] is a two dimensional crystal of carbon atoms, which form a hon- 
eycomb lattice with two distinct sublattices (A and B). The first Brillouin zone (BZ) 
has a hexagonal form (the blue area in Fig. lA), and the conduction band touches the 
valence band in six BZ corners [4] which form two non-equivalent triads of BZ corners, 
K and K' connected by the reciprocal lattice vectors, G and G'. Low-energy elec- 
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tronic excitations in the momentum space are located in the vicinities of the points 
K and K', i.e. belong to one of the two valleys with the linear 'Dirac-type' spec- 
trum, e{p) = ±vp where p is the momentum counted from one of the K-points and 
V ~ lO^cm/sec. Both the gapless spectrum and the valley degeneracy follow directly 
from the symmetries of the honeycomb lattice. In pristine graphene, the honeycomb 
lattice is stable against spontaneous structural changes. 

Recently, several types of adatoms were used to dope graphene in attempts to 
tailor properties of graphene-based devices [HI El [71 E] • Below, we consider theoreti- 
cally a particular example. Fig. 2A, of adatoms, such as alkali atoms [HITO], Ca or 
Al whose stable positions are above the centers of the hexagons. A single adatom of 
this type preserves rotational and reflection symmetries but breaks the translational 
symmetry of the lattice. Therefore, it can scatter electrons between valleys. The 
intervalley scattering generates the Friedel oscillations (FO) of the electron density, 
which amplitude rather slowly decays when the distance from the adatom |Tl]. The 
pattern and period of such FO are determined by the wave vectors 6K. = ±(K' — K) 
transferred upon scattering between the valleys. Due to a peculiar relation of K, K' 
and the reciprocal lattice vectors in graphene, K = |(G + G') and K' = — |(G + G'), 
the FO formed around an adatom have a structure of a charge density wave with the 
hexagonal lattice pattern and the unit cell extended over three unit cells of graphene. 
In graphene with zero carrier density (pe = 0), the amplitude of these superlattice 
oscillations decays as inverse cube of the distance between the adatoms. The oscil- 
lations of the electron density induced by one adatom affect other add atoms, thus 
leading to the pairwise interaction between adsorbents, which is sensitive to their 



3 



position in the superlattice. In Fig. 2 (A and B), we compare the potential landscape 
for a probe adatom created by several others due to their FO, to stress that its am- 
plitude is substantially enhanced by ordering. Note that the long-range RKKY-type 
[T2] interaction in a low-density 'gas' of adatoms as well as the ordering it promotes 
has little to do with those in dense aggregates of adsorbents. Indeed, the interaction 
at atomic distances is mediated by local lattice deformations - phonons. Such an 
interaction decays exponentially as a function of distance between the adatoms, as 
opposed to the power law 1/ B? decay of the electron- mediated coupling. 

Figure 2 illustrates an example of hidden structural ordering of adatoms sprinkled 
over graphene. Without mosaic coloring (or a superlattice mesh) it would be diffi- 
cult to distinguish the ordered configuration of adatoms (Fig. 2D) from a disordered 
state (Fig. 2C). With the help of colors, one can identify a triple-size unit cell of 
the superlattice, with three non-equivalent adatom positions (red, blue, and green) 
between six carbons, which resembles a hexagonal Kekule-type lattice [Ml [TO] . 
The FO-mediated interaction Vij between two adatoms i and j depends on whether 
they occupy equivalent (same color) positions on the superlattice, or not. This con- 
sideration maps the problem of hidden Kekule mosaic ordering of adsorbents onto 
the three-state Potts model [15] with a random in strength 'ferromagnetic' coupling 
of species. 

To estimate critical temperature of the hidden ordering, one has to evaluate the 
function Vij{Yi — vj), where and are adatoms locations. We use the technique 
developed for the studies dealing with disorder in graphene [16l [171 [H] and describe 
electrons by four-component Bloch function ip (taking into account its sublattice 
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composition and valley degeneracy) and the 4x4 Hamiltonian, 



H = v{p- a)U, + Xhva ^ (u; • U)5{r - r^). (1) 

i 

Here, Pauli matrices (Tx,y,z act on sublattice indices and ^x,y,z on the valley indices 
of ip [ini [T71 [11] . The first term in Eq. ([1]) is a familiar Hamiltonian for Dirac-type 
electrons in pristine graphene. The second accounts for the intervalley scattering of 
electrons by adatoms, with A being the dimensionless coupling constant (realistically, 
A < 1 since hv/a is of the order of the bandwidth in graphene). Unit two-component 
vector Uj = (cos ^^y^, sin ^^j^) specifies which of the three non-equivalent positions 
the i-th adatom occupies on the superlattice, with = —1, and —1 (red, blue, 
and green hexagons). 

Using thermodynamic perturbation theory and the standard RKKY approach [12] 
we express the interaction between adsorbents mediated by electrons as 

/oo 
rfr(ui ■ n)G{vi - Vj, r)(ui ■ n)G{rj - r^, -r), 
-oo 

^'^ 47r(t;2r2 + r2)3/2' 

where r is imaginary time and G{r, r) is the zero-temperature Green function of 
Dirac-like electrons. Strictly speaking the equation 2 is valid at T = 0. However one 
can use it as long as T < hvp^^^. The, trace (Tr) is taken over the sublattice and 
valley indices. The electron spin degeneracy is accounted for by the overall factor of 
2. The resulting electron-mediated 'ferromagnetic' interaction between adatoms at a 
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distance Irj — r J ^ a, 



u,; ■ Uj 2Tc(mi — rrii) 



= ~'^ |r,-r,|3p3/2 ' u,-u, =cos , (3) 



2 a 



has a long-range tail, oc |rj — The typical interaction energy scale, J is the 

interaction strength at a mean distance between the nearest neighbors, ~ p~^/^ (recall 
that p ^ a"^). 

To evaluate the critical temperature, T^, we modeled the ordering transition nu- 
merically. We used the cluster Monte Carlo algorithm [18| to compute statistical 
moments 

A^ui,...,ujv \^-^i=l N 

^2n = =^ _J_y-y.. ' 

\ p 2T ^ "^'J 

^Ui,...,U]v 

for 10 realizations of quenched Poissonic distributions of = 2 x 10^ adatoms. 
The ordering transition can be detected by a sudden rise of the order parameter 
M = \/M2, from M(T > Tc) = to M(T < Tc) = 1 accompanied by decrease of 
f] = M4/M2, from ri{T ^ Tc) = 2 (set by the central limit theorem for a large number 
of uncorrelated clusters) to ri{T <^ T^) = 1. Results of the numerical analysis are 
presented in Fig. 3. The transition temperature turned out to be 



T, ^ 8J ~ 4X\a^pf' — . (4) 

a 



For example, for A ~ 1, just 1% coverage of graphene by adatoms should generate Tc 
in the room temperature range. 
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Since the mobility of adatoms on graphene strongly depends on temperature, the 
higher the adsorbent density p, the higher is, and the quicker the self-organization 
should establish upon cooling. Note that the aggregation of adsorbents, such as 
discussed in Refs. [191 120] , would be a much slower process in a dilute system. At 
T < Tc the proposed partial ordering suppresses adatoms diffusion leading to a further 
slowdown of aggregation. 

The value of may also depend on the concentration of electrons (or holes) in 
graphene. Finite carrier density leads to the additional modulation of the FO, with 
the period twice as small as the electron Fermi wavelength ~ pe [Hj- For Pe ^ P 
these modulations would make the sign of the interaction between adatoms random 
and, thus eliminate the ordering. Therefore it seems to be possible to control hidden 
ordering of adsorbents electrically, by filling or depleting the flake with carriers - the 
method already in use to fine-tune the ferromagnetic transition temperature in thin 
films of dilute magnetic semiconductors [2T] . 

The self-organization of an apparently random ensemble of adatoms into a Kekule- 
type ordered state drastically changes electronic spectrum in graphene. Adatoms 
that preferentially occupy one of the three equivalent positions in the supercells over 
a length scale L ^ a/1/p can Bragg scatter electrons between the two valleys coher- 
ently. This implies the Brillouin zone folding in Fig. 1 (B and C): all of the points K 
and K' of the original BZ are projected onto the F-point of a smaller BZ correspond- 
ing to the superlattice with a triple unit cell. Simultaneously, a gap. A, opens in the 
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electronic spectrum 



eip) = ±^/{vpy + A\ A = AaV— . (5) 

To derive Eq. one can substitute the second term of the Hamiltonian in Eq. 
([1]) by its average, for example, for all adatoms positioned on yellow hexagons, and 
diagonalize the resulting matrix, H = vcr ■ p + XvapUx- 

One can think of several ways to experimentally detect the hidden Kekule mosaic 
order. One is to use the angle- resolved photoemission spectroscopy (ARPES). The 
latter technique is not only a natural method to reveal the formation of the spectral 
gap. It can also provide a direct evidence of the BZ folding. Indeed, ARPES measures 
simultaneously the energy and all three momentum components of the photo-emitted 
electrons. While at low energies only the vicinity of the BZ corners K and K' can be 
seen in pristine graphene [221 123] , the Bragg scattering by the self-organized adsor- 
bents generates an ARPES signal also at the F-point of the BZ in Fig. lA. Another 
signature of the Kekule ordering would be a bright appearance of the D peak in the 
phonon Raman scattering [23]: the excitation of a BZ phonon forbidden by momen- 
tum conservation in pristine graphene. Finally, the gap in the electronic spectrum 
would dramatically affect charge transport in graphene. This may offer numerous 
opportunities for graphene-based electronics. 
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Figure legends 



Figure 1: (A) Brillouin zone of graphene in the reciprocal space is shown as a blue 
hexagon. The two valleys K and K' are situated in the six corners of the hexagon, 
which are identified via reciprocal lattice translations generated by vectors G and 
G'. (B) Brillouin zone folding due to the ordering transition. The folding leads to 
the identification of the valley points K and K' with the F point in the center of 
the Brillouin zone. (C) The energy surface in the folded Brillouin zone. Due to the 
interaction between the valleys a gap opens in the spectrum. 



Figure 2: Kekule mosaic ordering of adatoms of the same chemical element on 
graphene lattice. Panels (A) and (B) show the potential landscape that an extra 
atom would see in the presence of four atoms already shown. Coloring of the atoms is 
introduced to reveal their position within the Kekule superlattice, as shown in panels 
(C) and (D). From a comparison of (A) and (B) one can see that adatoms placed on 
unicolor tiles enhance the potential landscape forcing other atoms to occupy tiles of 
the same color. 
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Figure 3: The Kekule mosaic order parameter M as a function of temperature. The 
phase transition to the ordered state is characterized by a rise of M accompanied 
by a sharp drop of the measure of finite-size fluctuations M4/MI. The data were ob- 
tained using the cluster Monte Carlo algorithm. The error bars indicate the standard 
deviation of the thermodynamic quantities in an ensemble of 10 random realizations 
of Poisson distributions of = 2 x 10^ atoms in the plane. 
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